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Abstract 

We study the well-posedness of the Hele-Shaw-Cahn-Hilliard system modeling binary 
fluid flow in porous media with arbitrary viscosity contrast but matched density between 
the components. For initial data in , s > f -1-1, the existence and uniqueness of solution 
in C([0, T]; iJ^) nL2(0, T; 7?^+^) that is global in time in the two dimensional case (d = 2) 
and local in time in the three dimensional case {d = 3) are established. Several blow- 
up criterions in the three dimensional case are provided as well. One of the tools that 
we utilized is the Littlewood-Paley theory in order to establish certain key commutator 
estimates. 

1 Introduction 

The modeling and analysis of multi-phase fluid flow is a fascinating, challenging and important 
problem [18^ [5]. Well-known two phase fluid examples include the coupled atmosphere-ocean 
dynamical system with water and air being the two phases, as well as the system describing 
displacement of oil by water in oil reservoir (usually porous media) [5j . 

A common approach to two phase flow that are macroscopically immiscible is the sharp 
interface approach where the two phases are separated by a sharp interface T(t). In the case 
of flow in porous media, the dynamics of the system is then governed by the two phase Hele- 
Shaw (Darcy) system (Muskat problem) [201 \T7\ [25] together with two interface boundary 
conditions: (1) continuity of the normal velocity; and (2) pressure jump proportional to the 
(mean) curvature. The normal velocity of the interface is set to be the normal velocity of the 
fluids. The local in time well-posedness of the sharp interface model with or without surface 
tension is known [SJ [3l [13]. Global in time well-posedness with surface tension [Ml [10] and 
2D without surface tension [27] is also known under the assumption that the initial data is 
a small perturbation of a flat interface or a sphere. Nevertheless, the sharp interface model 
encounters serious difficulty with physically important topological changes of the interface 
(possibly undefined curvature), especially in terms of pinchoff and reconnection that are 
important in applications [ij [20]. 
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As an alternative approach, one could consider the so-called phase field models (or diffuse 
interface models) where an order parameter c is introduced and a capillary stress tensor is 
used to model the interface between the two fluids and the forces associated [1]. The sharp 
interface is then replaced by a thin transition layer and hence we avoid the difficulty of 
discontinuity. In this paper, we will consider phase field approach to two phase fluid flow 
with matched density in a Hele-Shaw cell or porous media. The dynamical equations are 
given by the following Hele-Shaw-Cahn-Hilliard system |20^ I12j: 

' V-M = 0, 

ct + u-Vc = ^An, 
, c(0,x) = co{x), 

where u is the fluid velocity, c is the order parameter which is related to the concentration 
of the fluid, the chemical potential fi depends on the order parameter c and is given by 

Mc) = /o(c) - CAc, (1.2) 

and Pe is the diffusion Peclet number, C is the Cahn number, and M is a Mach number. 
Furthermore, r](c) is the kinematic viscosity coefficient satisfying 

r7GC°°(R^), < A < 7?(c) < A < oo, (1.3) 

the Helmholtz free energy /o(c) is given by the classical double well potential 

/o(c) = (c^ - If. (1.4) 

In the above system (II. ip . p is not the physical pressure but the combination of certain 
generalized Gibbs free energy and the gravitational potential (see ^20j for more details). 
This model can be also viewed as the Boussinesq approximation of more general model with 
arbitrary viscosity and density contrast [20]. One may formally recover the sharp interface 
model by taking appropriate limit within the Hele-Shaw-Cahn-Hilliard system (jl.ip [20]. 
We will assume that the fluid occupies the two or three dimensional torus T"', d = 2,3 for 
simplicity. 

Besides applications in two phase flow in porous media and Hele-Shaw cell, certain simpli- 
fled versions of this HSCH model has been also used in tumor growth study [30]. Moreover, 
unconditionally stable schemes has been developed ^29j and the existence of certain type of 
weak solutions (without uniqueness) is also derived [15] for the case with matched density 
and viscosity. 

The goal of this manuscript is to study the well-posedness of the matched density Hele- 
Shaw-Cahn-Hilliard system (jl.ip with arbitrary viscosity contrast. 

The Hele-Shaw-Cahn-Hilliard system can be formally viewed as an appropriate limit of 
the classical Navier-Stokes-Cahn-Hilliard system [H [20l [T6] which is a popular phase field 
model for two phase flow although no rigorous justiflcation is known yet. There are a lot of 
works on the Navier-Stokes-Cahn-Hilliard system including local in time well-posedness in 2 
and 3 dimensional and global in time well-posedness in 2D under various assumptions [U [7] . 
In fact the global in time well-posedness of the 2D Navier-Stokes-Cahn-Hilliard system is 
recently resolved [l] using a very different set of tools than employed here. Mathematically 
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speaking, the difficulty associated with the Hele-Shaw-Cahn-Hilliard is about the same as 
those associated with the Navier-Stokes-Cahn-Hihiard: we gain the advantage of dropping 
the nonhnear advection term in the velocity equation but also lose the regularizing viscosity 
term; and their scaling behaviors are very similar. We refer to |21[ l22l l23l S] and references 
therein for more related works on the Navier-Stokes-Cahn-Hilliard system.. 

The rest of the paper is organized as follows. We prove a key estimate on the "pressure" 
in the second section. This estimate is nontrivial due to the variable coefficient introduced 
with the mismatched viscosity. New estimates on certain commutator operators in fractional 
derivative spaces are needed and they are derived in the Appendix. In section three we 
present the local in time well-posedness based on certain modified Galerkin approximation 
of the HSCH system and the "pressure" estimate from section 2. In section 4 we provide a 
Beale-Kao-Majda type blow-up criterion and prove that the system is global in time well- 
posed in the two dimensional case. We provide a refined blow-up criterion in the 3D case in 
section 5. 

2 The estimate of the pressure 

In this section, we present the estimate of the modified pressure p. Taking the divergence for 
the second equation of (II. ip . we find that 

divf^rVp) = divf^-T/^WVc) =^divF(c). (2.1) 
\?7(c) / ^^?(c) ^ 

This variable coefficient problem is dealt with utilizing commutator estimates that we derived 
in the Appendix. The commutator estimates themselves are derived utilizing Littlewood- 
Paley decomposition. 

Proposition 2.1 Let s > and c € H^'^'^{T'^), and p be a smooth solution of ^2.1\) . Then 
the solution p satisfies 

llVplliT. < ^(||c||loo)(1 + ||Vc|U-) (1 + \\c\\H2f\\c\\Hs+2. (2.2) 

Here k = [2s] + 1 and J- is an increasing function on R"*". 

Proof. Thanks to (jl.3p . a straightforward energy estimate yields that 

\\Vp\\l2 < C||/i(c)||i2||Vc||Loo < C7(l + ||c||ioo)||Vc||Loo||c||j^2. (2.3) 



Taking the operator {D)^ to (j2.ip to obtain 




= dw{A + B), 



from which and the energy estimate, we infer that 

WVpWhs <C{\\A\\L2 + \\Bh2). 
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Due to the definition of /u(c), we have 

^ /i(c)Vc = ^/^(c)Vc - C^AcVc = V(7i(c) - AcVg2{c), 



ri{c) ri{c) ri{c 

for some 51,52 with gi(0) = §2(0) = 0. We have by Lemma [631 that 

\\{DrVgi{c)\\L2 < J'i\\c\\L^)\\c\\Hs+., 

and using Bony's decomposition to write 

AcV52(c) = TAcVg2{c)+R{Ac,Vg2{c)) 

= divrvcV(72(c) - Tvc • VV52(c) + i?(Ac, V52(c)), 

then from the proof of Lemma 16.21 it is easy to see that 

i|(L>)^AcV<72(c)||L2 < Ti\\c\\L^)\\Vc\\Lo.\\c\\Hs+2. 

Thus we obtain 

<-^(||c||L-)(l + ||Vcb^)||c||^,.+2. 
and by Lemma [631ES and for s G (0, 1], 

\\B\\l2 < •7^(I|c||l°°)I|c||h»+2||V]5||l2 < J"(||c||Loo)||Vc||Loc||c||^^2||c||jys+2. 

Thus we obtain that for s € (0, 1], 

IIVpIIj^s <-F(||c1|l^)(1 + ||Vc||l-)(1 + ||c||j^2)|1c||h.+2, (2.4) 



For general s, we will prove it by the induction argument. Let us assume that for s € 
2 ' 21 



|], we have 



||Vp||//s < J-(||c||loo)(1 + ||Vc||l^)(1 + ||c||^^2)''||c||^^.+2. 

Note that (12. 4p means that the cases of A: = 1,2 hold. Now let us assume s G (|, ^y^]- We 
infer from Lemma 16.41 and Lemma 16.31 that 

\\B\\l2 < Ti\\c\\L^){\\c\\Hs+2\\Vp\\L2 + ||c||j,2||Vp||^._l). 

Then from ()2.3I) and the induction assumption, it follows that 

\\B\\l2 < -F(||c||l-)(1 + ||Vc||l.o)(i + ||c||^^2)^'+^||c||hs+2. 
Thus for s G (|, we have 

W^pWh' <7-(||c||loo)(i + ||Vc||loo)(i + |1c|1j^2)^+'||c||h.+2. 
This completes the proof of Lemma 12. li ■ 
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Remark: Instead of relying on the estimates from the appendix which depend on the Littlewood- 
Paley theory, classical energy method might work as well if we are content with less sharp 
and less general results. For instance, if Vc G L°°{T'^) and c € H^{T^) for k G Z+, classical 
elliptic estimates may lead to 

\\Vp\\h^<C\\F{c)\\h^, 
where C depends on ||Vc||loo and ||c||j|^fc. And a straightforward product estimate gives 

It-F'Cc)!!//;. < C'(||c||L°°)(||c||j:/fc+l + |[Vc|[l«= ||c||j:^fc+2 + ||Ac||L°°||c||j:^fc+l). 

This estimate is enough to obtain the local well-posedness of the system (jl.ip and global 
well-posedness in the 2D case in the space of 

c € C([0, T]-H^{T'^)) n L2(o, T; H^{T'^)), u G C7([0, T\-L^{T'^)) n ^^(o, T; ^^(T'^)) 

when combined with the L°°{H'^) fi L?'{H'^) a priori estimates from (Theorem 14. ip for initial 
data in , k > 2. However, in order to obtain the sharp blow-up criterion which in partic- 
ular implies the global existence of the 2-D system in general Sobolev spaces as specified in 
Theorem 13. H we need to establish the refined pressure estimate (12. 2p . Notice that (12. 2p is 
established for general (Hilbert) Sobolev spaces, and only a linear in ||Vc||ioo factor in the 
estimate appears in contrast to pure energy estimates. 



3 Local well-posedness 

In this section we prove the local well-posedness of the Hele-Shaw-Cahn-Hilliard system. The 
procedure is mostly standard except the pressure estimate. 

Theorem 3.1 Let co(x) G H''{T'^) for s > f + 1. Then there exists T > such that the 
system has a unique solution (c, u) in [0, T] with 

c e C{[0,T]; H' {T'^)) n L'^{0,T; H'+^iT'^)), u G C{[0,T]; H'^^iT"^)) n L'^{0,T- {T'^))- 
and satisfying the following energy estimate 

Mt)\\h+ [ \\c{T)\\l^+2dT < \\co\\h^ exp { [ G{T)dT). (3.1) 
JO Jo 

for t G [0, r], where 

G{t) = H\\c\\l^){1 + \\Vc\\lo^)\\\Vc\\l^ + ||c||^)'(l + ||c||h2)'(['^]+'\ 

Proof. We will use the energy method to prove Theorem 13.11 
Step 1. Construction of an approximate solution sequence. 

The construction of the approximate solutions is based on Galerkin method. Let us define 
the operator P„ by 

Pnfix) = Ae^-'^-, f,= [ /(x)e-2-'=-^dx. 

\k\<n 
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(3.2) 



Then we consider the foUowing approximate system of 

V • = 0, 

^ ~12ri[PnC„) {^Pn - PnCn) , 

dtCn + Pn{Un ' VP„C„) = ^ AP„;u(P„C„) , 
Cn(0,x) = PnCo(x). 

It is easy to see that 

||AP„/i(P„4) - ^Pn^l{Pncl)\\L^ < C{n, 11411^2, ||c2 ||i2)||c^ - cI\\l2. 

Taking the divergence to the second equation in ()3.2p gives 

Thanks to (jl.Sp . straightforward energy estimate yields that 

||Vp„||i2 < C{n, ||c„||j^2)||c„||i2, 
thus we infer from the second equation of (|3.2p that 

||^in||L2 < C{n,\\cn\\L^)\\cn\\L2■ 
Theve^oTe, we have 

\\Pn{ui ■ VPnci) - Pn{ul ■ VP„c2)||^2 < C(n, \\cI\\l2, || ^2 ) - cI\\l2. 

Thus, the Cauchy-Lipschtiz theorem ensures that there exists > such that the approx- 
imate system (I3.2p has a unique solution Cn G C([0, T„]; L^(T'^)). Note that P^ = Pn, PnCn 
is also a solution of (|3.2p . So the uniqueness implies that PnCn = Cn- Thus, the approximate 
system ()3.2p reduces to 

r V-Un = 0, 

^« = -T2i^i'^Pn - M/^(Cn)Vc„), ^^^^^ 
5tC„ + Pn{Un ■ Vc„) = ^AP„/i(c„), 

^ c„(0,x) = P„co(a;). 

In what follows, we denote T* by the maximal existence time of the solution Cn- Due to 
Pn Cn — Cn, the solution Cn is in fact smooth. 

Step 2. Energy estimates. 

Although the HSCH model (jl.ip has a natural energy (which is somewhat equivalent to 
estimate, see [20^ [29] and section 4 below), it is not sufficient for the strong solution. 

Therefore we have to derive estimates in Sobolev spaces with higher derivatives. 

For this purpose we take the H'^{T^) inner product of the third equation (13. 3p with c„ 

and obtain 

^ ^ II Cn IIh^ - ^(AP,t/i(c„),Cn)^, = -{un-VCn,Cn)jjs- (3.4) 
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Due to (|l.'2p . we see that 

-{APnlJ,{Cn),Cn)jj, = C||Ac„||^. - (A/q(c„) , C„) . 

We deduce, thanks to Lemma 16.21 that 

Cn ||l°° ) II*^™ \ \h^ 

||Ac„||//s. (3.5) 

and by Lemma 16.21 with a = 1, 

\{un-VCn,Cn)fjs\ < H^^n • Vc^ ||c„ 

< C[\\Un\\H=\\^Cn\\L-^ + \\Un\\^d_^\\VCn\\H''+i)\\Cn\\H^- (3.6) 

Thanks to (|3.3p . we find that 

||n„||//= < C'dl— ^Vp||h= + ||— ^//(c„)Vcn||//»). (3.7) 

By Lemma [6.2( Lemma 16.31 and Proposition 12. H the first term on the right hand side of (13. TP 
is bounded by 

-F(||c„||L-)(||c„||^.+ .||Vp||i2 + llVpll^^.) 

< •^(l|Cn||L°°)(l + ||Vc„||l°°)(1 + ||c„||j:^2)'^'*''^^||c„||j:^s+2, 

and the second term is bounded by 

-^(||Cn||L°°)(l + ||Vc„||loo)||c„||j:^s+2. 

Thus we obtain 

Iknlk'' < -^(|ICn||L°°)(l + |IVc„||loo) (l + ||c„ ||j:^2) '^"^'^^ ||c„ ||j:/s+2 , 

and similarly, 

\\Un\\^d_-^ < J'(||c„||l°o)(1 + ||Vc„||loo) (l + ||c„||^2)'^ ^\\Cn\\^d^i, 

from which and ()3.6p . we infer that 

\{Un ■VCn,Cn) < -F( || C„ || ^oo ) (l + || Vc„ || L-^ ) 

x(||Vc„||loo + ||c„||j^| )(l + ||c„||j:^2) " ||c„||j|^s+2||Cn||H«- (3.8) 

Here we used the following interpolation inequality: 

2 d d -j^ 

||Cn||^d^l ^ ||Cn II j:^2^ ||Cn |||f3 • 

Plugging (I33|) and I^Mj into ([331) yields that 
1 d 2 C 2 



< 7'(||c„||l=o)(i + ||Vc„||loo) (||Vc„||loo + ||c„||^| ) (l + ||c„||//2)'^'*^^^||c„||i^s+2||c„||H 



which along with Young's inequaUty imphes that 



d II ||2 , II ||2 

|2 



< J'(||c„||l«>)(1 + ||Vc„||lcx>) (||Vc„||l°o + ||c„||^?, ) (l + ||cn||H2) ||c„||h=. 
Then Gronwah's inequaUty apphed gives 

E'^it) =^ ||c„(t)|||,. + f \\Cn{T)\\jj.+2dT < \\co\\h^ exp ( f Gn{T)dT) (3.9) 

Jo Jo 

for t G [0, T^), where 

Gn{t) = ^(||c„||l-)(1 + ||Vc,,|Uco)'(||Vc„||loo + ||c„||^)'(l + ||c„||h2)'^''''+'^ 

Step 3. Uniform estimates and existence of the solution. 
Let us define 

f: t%up{t G [0,T:) : i?^(r) < 2||co|||^. for r G [0,t]}. 
From (j3.9p and Sobolev embedding, we find that 

Em < ||co||j/=exp(^(||co||i^O /V + llcWII^-a'Mr) 

Jo 

< ||co||//.exp(^(||co||iy.)(* + i^)), t^%T*). 
Here A{-) is some increasing function. Take T be small enough such that 

exp {A{\\co\\h^){T + T-2)) < I 

Now we can conclude that T* >T. Otherwise, we have 

K{t)<l\\cofHs for te[o,f:], 

which contradicts with the definition of T*. Thus the approximate solution (c^, itn) exists on 
[0, T] and satisfies the following uniform estimate 

\\Cn{t)\\h+ f l|Cn(r)||^.+2dT<2||co||H. (3.10) 

Jo 



for t G [0,T]. On the other hand, it is easy to verify from the third equation of (|3.3p 
that dtCn is uniformly bounded in L^{0,T; H^~^{'T'^)). Thus, Lions-Aubin's compactness 
theorem ensures that there exist a subsequence {cn,^,Unf.)k of (c„,u„)„ and a function c G 
L'^{0,T;H'{T'^))nL^{0,T;H'+^{T'^)) and w G L'^{0,T; H'-^iT^)) n L'^{0,T; H'iT^)) such 
that 

Ts' + 2/rpd\ 



in L'{0,T;H'+'{T'')), 



Un.^u, in L\0,T-H'' {T'')), 
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as k ^ +00, for any s' < s. Then passing to limit in (j3.3|) . it is easy to see that (c, u) satisfies 
(jl.ip in the weak sense and (c, n) satisfies ()3.ip . 

Step 4. Continuity in time of the solution. 

Revisiting the proof of ()3.9p . we can in fact obtain better estimate for c„ (thus for c): 

ll'^ll|°°(0,T;i?''(T'')) = X] 2^^1l^i^!li°°(0,T;L2) ^ C*' 

i>-i 

which will imply c G C([0, T]; ff'^(T'^)). In fact, for any e > 0, take N big enough such that 

2^''*|[AjC|||oo(o,T;L2) ^ 7- 

For any t € (0, T) and (5 such that t + (5 € [0, T], we have 

AT 

|jc(t + 5) -c(t)||^. < ^ 22i«||A,c(t + (5) - A,-c(t)||i2 + 

i=-i 

AT 



e 



Thus for I (5 1 small enough, we have 

\\c{t + 5)-c{t)\\'^ <e. 

That is, c(t) is continuous in H'^{T'^) at the time t, thus so does u. 
Step 5. Uniqueness of the solution 

Assume that (ci, ni) and (c2, 1*2) are two solutions of (jl.ip with the same initial data. We 
introduce the difference of two solutions: 

6c = Ci-C2, 5u = Ui-U2. 

Then {5c-,5u) satisfies 

dt5c + ui ■ V(5c + 5u ■ Vc2 = ^A(/i(ci) - 11(02)), 

= tt)tcl] (^Pl - M/^(^l)Vci) - U^^i'^iPl -P2) - I^(/X(C1)VC1 -^(C2)VC2)), 

^ SciO) = 0. 

Making L^(T'^) energy estimate yields that 

< C(||A(5c||l2 + ||(5„||i2)||(5c||L2. 
On the other hand, we can deduce from the equation of 6^ that 

115^.11^2 < C7(||(5,||i2 + ||V(pi -P2)||l2 + 
< C{\\5c\\l^ + \\A5c\\l^). 
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Thus we obtain 

j^\Mh<c\Mh, Pc(o)||=o, 

which along with Gronwall's inequahty imphes 5c = 0, and the uniqueness follows. I 

4 Blow-up criterion and global existence in 2D 

In this section we prove a Beale-Kato-Majda type blow-up criterion [24J for the Hele-Shaw- 
Cahn-Hilliard system. As an application, we obtain the global well-posedness in 2D. 

Theorem 4.1 Let co(x) € H'^{T'^) for s > | + 1, and {c,u) be a solution of stated in 

Theorem \3.1[ Let T* be the maximal existence time of the solution. IfT*< +oo, then 



I |lVc(t)||ioo(it = +00. (4.11) 
Jo 



In particular, this implies T* = +oo for d = 2. That is, the system is globally well-posed 
in 2D. 

Proof. First of all, we derive the basic energy law of the system. Multiplying by ^ on both 
sides of the third equation of (jl.ip . we get by integration by parts that 

I ctixdx + I u ■ V cjidx = — —- I \Vfi\'^dx. 

Due to the definition of /i, we have 

/ ctndx = ^([ fo{c)dx + ^ [ iVcpdx), 

J'J^d at \ J'J^d Z J ^d / 

and due to V • n = 0, 

/ w ■ Vc/xdx = -M / u ■ {Vp - ^fiVc)dx = 12M [ ri{c)\u\'^dx. 

Jrpd Jrj^d M Jrpd 



Thus we obtain the following classical energy equality |20j 

/ fo{c)dx + — / iVcpdx) + ^ f \V^JL\^dx + 12M / r]{c)\u\'^dx = 0. 
. J '^d 2 J'j^^ ^ Pe J'j^d J 



That is. 



where 



E{t) + ^ f f \Vfi{T)\'^dxdT + 12M [ [ r]{c)\u{T)\'^dxdT = E{0), (4.12) 
E{t)=^[ fo{c{t,x))dx + ^ [ \Vc{t,x)\^dx. 
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From the energy equality (|4.12p . it follows that 

imwii + ^ 1^ wvfiwi^dT < E{o). 

On the other hand, we have 

||VAc||i2 < C{\\S7fih2 + ||Vc|li2 + \\c^Vc\\l2). 
and by Sobolev inequality, 

\\c'Vc\\l2 < CllclllellVcll^e <C|Ic||^i||c||;,2 

< C||c|||,||c|||3 <C||c|||,i + ^||c||h3, 

which implies that 

< C(I|V/x||l2 + ||c||^^i + ||c|| 

Therefore we conclude that 

I|c||l-'{0,T;H1) + I|c||l2(0,T;H3) < C{T, ||co||hi)- (4.13) 

Next, we derive energy estimate of the solution. We have 

^1"^^"^^ + ^"^^'^"^^ = -(^•Vc,A2c)^, + ^(A/^(c),A2c) 

< \\u\\l2\\Vc\\l^\\A''c\\l2 + ^||A/^(c)||i2||A2c||i2.(4.14) 

It is easy to verify that 

< C{\\Vp\\l2 + Mc)Vc\\l2) 

< C(||Vc||L^||Ac||i2 + (||c||i3 + ||c||i.)l|Vc||i6) 

< C{\\Vc\\l^ + |lc|li,3 + llcllls) \\c\\h2, 

and 

l|A/^(c)||z.2 <C7(l + ||c||i^)||c||^2. 
Plugging them into (j4.14p yields that 

|||Ac||i. + IIA^dli, < C{1 + llVclli. + llclli. + llcllia + \\c\\'3)\\c\\jj., 
which along with Gronwall's inequality leads to 

||c||h2 < ||co|Ih2 exp (C / H{T)dT), (4.15) 

where H{t) = 1 + \\Vc\\l^ + ||c||i^ + ||c|||3 + ||c||ii 
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Now we are in position to prove the blow-up criterion. We will prove it by way of 
contradiction argument. Assume that T* < +00 and 

/ ||Vc(t)||ioodt < +00, 
which together with (j4.13p and Sobolev's inequality implies that 



I H{T)dT < +00, 
^0 



for example, 



'0 ^0 
Then we infer from (j4.15p that 



c{t)\\fgdt <C \\c{t)\\]^i\\c{t)\\H3dt < +00. 



\l°°{Q,T*;H'^) < +00, 



which implies that 



/ G{t)dt < +00, G{t) be as in Theorem 13.11 
Jo 



Then the energy inequality (jS.ip ensures that 

\m\\h+ r 

te[o,r 



sup ||c(t)|||^a + / ||c(T)||^s+2dr < +00, 
:fo,r*i JQ 



which means that the solution can be continued after t = T* , and thus contradicts with the 
definition of T*. 

As an application of blow-up criterion, we can deduce the global existence in 2D. Indeed, 
in two dimensional case, we get by Gagliardo-Nirenberg inequality and (j4.13p that 



[ \\Vc{t)\\Udt<C [ ||c(t)||^i||c(t)||^3rft < +00, 
Jo Jo 

which implies T* = +00 by the blow-up criterion. 

5 A refined blow-up criterion in 3D 

We first turn to a simple model relating to the Hele-Shaw-Cahn-Hilliard system: 

( u = - Vp + AcVc, V • li = 0, 
\ Q + n • Vc + A^c = 0. 

For this system, we still have the energy equality: 



(5.16) 



Vcit)\\i,+2 / ||VAc(T)||i, + ||n(r)||i,dr = ||Vco|| 



L2. 
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Moreover, if c is a solution of ()5.16p . then c\{t,x) c(A*^t, Ax) is also a solution. It is easy 
to see that 

\\Vcxit,x)\\L2=X-2-^\\S/c{XH,x)\\L2, / \\VAcx{T)\\l2dT = X^-^ \\VAc{T)\\l2dT. 

Jo Jo 

Thus, the energy is scaling invariance for d = 2. From this view of point, the 2D system is 
critical and the 3D system is supercritical like the 3D Navier-Stokes equations. Due to the 
bi-Laplacian A^, there is no maximum principle for this system, which is the main obstacle 
to obtain the global existence in 3D case. For the 2D critical QG equation 

et + {-A)'^e + u-ve = {), u= {- {-A)-'2d^,e,{-Ay^d,,e), 

Caffarelli and Vasseur [5] proved the global regularity of weak solution. The key step of 
their proof is to prove the Holder continuity of the solution by using the DeGiorgi method. 
Note that the QG equation has maximum principle. For the 3D Hele-Shaw-Cahn-Hilliard 
system, we also show that the Holder continuity of the solution will control the blow-up of 
the solution. 

Theorem 5.1 Let a G (0, 1) and co{x) G F''(T^) for s>3. A ssume that (c, u) he a solution 
of (E2P stated in Theorem \3.1\ Let T* be the maximal existence time of the solution. Lf 
T* < +00, then 

j-T* s 

/ \\c{t)\\gcdt = +00. 
Jo 



Proof. We will prove it by contradiction argument. Assume that T* < +00 and 

i^^dt < +00. (5.17) 



Jo 
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Taking Aj to the third equation of (|1.1|) to obtain 



dtAjC + p^A^AjC = -Aj{u ■ Vc) + — AAj/^(c) 



Making L (T ) energy estimate, we get by Lemma l6. II that for j > 0, 

|||A,c||i. +c24^||A,-c||i. < C{\\Aj{u.Vc)h2 + || A/^(c)||,.2) || A,c||^2. 
Dividing the above inequality by ||Ajc||^2 gives 

^||A,-c||,.2 + c2'^\Ajch2 < C{\\Aj{u ■ Vc)\\l2 + ||A/^(c)||,.2), 
which implies that 

\\Ajc{t)U2 < ||A,co||z.2 +C fe-''''^'-^\\\A,{u-Vc){T)h2 + ||A/^(c(r))||i2)dT. (5.18) 

Jo 
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We denote 

Mb^^ =^ sup 2^'\\Ajc\\L2. 

i>-i 

Using the definition of Sobolev space, it is easy to find that 

||c||^._. < 2-2-^||c|||.^ < C\\c\\ls^, Ve > 0. 
i>-i 

It fohows from ([5J8]) that 

||c(t)||B3^ < \\c{t)\\L2 + WcoWh^ 

+C7sup23^- re--24^(*--)(||A,(n.Vc)(r)||i. + ||A/^(c(r))||^2)dr. (5.19) 
i>o Jo 

Now we claim that 

||Aj(n- Vc)||i2 < C2J'(^-")|| u||i2||c||c"^. (5.20) 

Now we have 

||n||i2 < C\\h{c)Vc\\l2 < C||c||i^3-c«||c||c- +C(||c||i3 + ||c||^6||c||loo)||Vc||z,6 
< C{'^ + \\4m + I|c||^i)||c||c"||c||s3^. 

Here we used the product estimate 

||AcVc||i2 < C||c||^^3-c||c||c- < C\\c\\b3 ||c||c", 

z.oo 

which can be proved as in Lemma 16.21 And similarly we have 

llA/^(c)||i2 <C7(l + ||c||^.)||c||^2. 
Plugging the above estimates into (j5.19p yields that 

umBi < iic(t)iiL2 + iico||h3 

+Csup2^-(4-") re-^2^^(*-)(l + ||c||^i + ||c||^i)(l + ||c||2,.)||c||o3 dr, 
j>0 Jo •°° 

which along with Holder inequality gives 

l|c(i)|lL°°(0,t;B3^) < \\c{t)\\L°°{0,t;L2) + I|co||h3 

+ (1 + l|c||L-(0,i;Hi) + \\4l^(0,t;H^)) (*~ + 11^11^1 (o,4;C-)^ (0,t;B3 ^) • 

The above argument is still valid on the interval [T,T*) for T < T* . Thus we get by using 
(ITOjl that 

lic(OllL°°(T,T';B|^) < IWoWm + \\co{T)\\h3 

+C(||co||HO((r*-T)U||c||^^^^^^,^^^pi|c||,.(^,^.^3^). 
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Due to (|5.17p . we can choose T such that 

C{\\co\\m){iT* -T)^ + \\c\W ) < I, 

Then we obtain 

I|c(OIIl-(t,t*;B|^) < 2(||co||hi + ||co(r)||H3), 
which imphes by ||Vc||ioo < C||c|| ns that 

2,00 

/ ||Vc(t)||^oo(it < +00, 

which is impossible by Theorem 14.11 if T* < +00. 

It remains to prove (|5.20p . As in proof of Lemma 16.21 we have 

Aj(n-Vc) = Aj ^ • VAfcC + Aj ^ A^n • VS'fc^ic 

|j-/c|<4 |j-fc|<4 

+Aj ^ AfcU- VAfc'C = ^l+^2 + ^3• 

fc-fc'|<l,fc>j-3 

We get by Lemma [6T] that 

Pi||l2<C ^ ||5fc_in||i2||VAfcc||L.o <C2^-{i-")|| 

^IIl2||c||c", 

|j~fc|<4 

and for ^42, 

P2IIL2 < C l|Afcu||i2||V5fc_ic|U«. 

|i-fc|<4 

< C||^x||^2 E 2'||A,c||i^ 

|j-fc|<4f<fc-2 

|j-.fc|<4^<fc-2 

and due to V • n = 0, 

11^3 < II Yl V- (AfcuAfc,c)||i2 

|A:-fc'|<l,A:>j-3 

< C2^ Y 2-'='"||n||i22'='"||Afe/c|[L^ 

\k-k'\<l,k>j-3 

< C72J'(i-°)||n||i2||c||c-. 

Then the inequahty (I5.20p fohows from the estimates of Ai,A2 and ^3. The proof of Theorem 
15. H is completed. ■ 
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6 Appendix 



Let us first recall some basic facts about the Littlewood-Paley theory. Let 99, x be two 

:i<ici<|}, suppxc{iei<|; 



functions in C°°(T'^) such that supp^ C{|<|,^|<|}, suppx C {|^| < 1} and 



X(0 + E^(2^'^) = l- 

J>0 



Then the Littlewood-Paley operators are defined by 



/ T 

i-i 



fc=-i 

Some classical spaces can be characterized in terms of Aj. Let s € R, the Sobolev space 
H'{T'^) is defined by 

H^T") {u G V'{T^) : ||n||l,. ^ 2'n\^ju\\l, < 00}. 

i>-i 

We denote by {u,v)h^ the inner product in H^{T'^). And for s € (0,1), the Holder space 
C'{T'^) is defined by 

C'{T^) 'M{u e V'{T'^) : ll-ullc^ =^ sup 2^'\\Aju\\l^} . 

i>-i 

We refer to [28] for more details. Let us recall Bony's decomposition from [6]: 

fg = Tfg + TJ + R{f,g), (6.1) 

where 

Tf9=Y.^^-^f^i9, R{f,9)= ^jf^ra- 

i>-i li-i'l<i 
We also denote R{f,g) = Tgf + R{f,g). 

Lemma 6.1 IffJ Let /c€N, 1 < p < q < 00. Then there exists a positive constant C 
independent of j such that 



|A,/||iP < C2-^'^ sup ||9"A,/||lp, j > 0. 

\a\=k 



Lemma 6.2 Lei s > 0. Then there holds 

\\fg\\Hs < C{\\f\\L^\\g\\H^ + \\f\\H^\\g\\L^). (6.2) 
IfO<a<^, then there holds 

WfgWHs < C{\\f\\Hs\\g\\L^ + ll/ll^|_J|ff||H^+0- (6-3) 
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Proof. The inequality (j6.2|) is classical, see [T^. Here we only present the proof of (|6.3p . 
Using the Bony's decomposition (j6.ip to write 



A,{fg) = AjiTjg) + A,{Tgf) + A,Rif,g). 
Taking into considering the support of Fourier transform of the term Tfg, we have 

Ii'-il<4 

Due to < o" < i, this gives by Lemma l6. II that 



which implies that 



fc<i-i 



l^j(^/fi')|lL2 ^ X] ll'S'j'-l/llL°°ll^j'fi'|lL2 

li'-il<4 



li'-il<4 



(6.4) 



(6.5) 



< C2-^%-||/||^._J|<7lb=+^, 

here and hereafter {cj}j>-i denotes a sequence satisfying ||{cj}j>_i||£2 < 1. 
Similarly, we have 

\\A,{Tgf)\\L2 < C Y \\Sj'-ig\\L^\\Aj,f\\L2 
li'-il<4 

li'-il<4 

< C2-^%-|b||L-||/k- 

Noticing that, after taking into account the support of the Fourier transforms, 

AjR{f,g)= Y A,{AyfAjng), 

j',j">j-3;\j'-j"\<l 

it follows from Lemma |6. II that 

\\A,R{f,g)h2 < C Yl 2^h\^j'f\\L4^j"9\\L^ 

j',j">j-3;\j'-j"\<l 

< C2-^^ Y 2(^--^-')(f+^)2^-'(f-'^)||A,./||^22^-"(^+'^)||A,-.5llL2 



i',i">i-3;|j'-i"|<i 



Thanks to the definition of Sobolev space, (j6.3p follows from (|6.4p - (l6.6p . 



(6.6) 
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Lemma 6.3 128'j Let s > 0. Assume that F{-) is a smooth function on R with F(0) = 0. 
Then we have 



||F(/)||H»<C7(l + ||/||i^)W+i 



where the constant C depends on sup ||F^^^(t)||ioo . 

fc<W+2,|t|<||/||L°° 

Lemma 6.4 Let s > 0. Then there holds 

WiOrifg) - fiDygh^ < C{\\f\\Hs^49\\L^ + ll/k^ ll^ll^.- 1 ) • 
If s & (0,1], then we have 

mrifg) - fiorgh^ <c\\f\\Hs+4gh2. 

Here the Fourier multiplier (D)^ is defined by 



Proof. Using Bony's decomposition (j6.ip to write 

{DYifg) = {DYiTfg) + {DyTJ + {DYR{f,g), 
f{DYg = TfiDYg + T^D)^gf + R{f, (DYg). 

Thus we have 

(DYifg) - f{DYg = {DYiTfg) - Tj{DYg + 7r(/,5), 

where 

7r(/, 5) = {DYTgf + {DYR{f,g) - T^DYgf - R{f, (DYg). 
As in the proof of (|6.3p , we can deduce by Lemma 16.11 that 

\Hf,g)\\L^<c\\f\\H^+2\\g\\L2. 

We illustrate the process by working out the estimate on the first term. Thanks to Lemma 
16. H we have 

UDYTJWh = W^jiDYTJWh <cY, 22^-^||A,-r,/||i. 

i>-i j>-i 

< C 2'^'\\Sy^ig^yf\\l, 

b-j'l<4 

< c 2''n\Sr-ig\\l^\\^rf\\h 

li-i'l<4 

< c Y 2^^-(^+i)|biii.iiA,./iii. 

Ii-i'l<4 

< C\\g\\l4f\\'.<C\\g\\l4f\\ Hs+2- 
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Let m{^i,^2) be the symbol of the paraproduct operator Tfg. Then {D)^{Tfg) — Tf{D)^g 
has the symbol 



m 



(ei,6)((6 + e2)'^-(6)0, 



which is supported in the region |^i + ^2] ~ I ^2 1- By the fundamental theorem of calculus we 
have 

Jo 

It is easy to verify that 6)V/i*(t6+6)(6)^"^"'* with £ [0, 1] is a Coifman-Meyer 

paraproduct uniformly for t G [0, 1]. Then we have 

\\{Dr{Tfg)-Tf{DrgU2 < C\\{D)'-' fUAiDr^'-'gh'^ 

for 61 G [0, 1], i + i = i and 1< g < 00, see P. 106 in |31j. Taking 9 = ^ q) = (00, 2) for 
d = 2, and 9 = 0, {p, q) = (6, 3) for d = 3, we obtain 



\\{DnTfg) -Tf{Dyg\\L2 <C\\f\\H49\\^s-^,- 
In case of s G (0, 1], taking 9 = 1 — s and {p, q) = (00, 2) to obtain 

\\{DY{Tjg) - Tj{Drg\\L2 < C\\f\\Hs+2\\gh2. 
This completes the proof of Lemma [61 
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